Analysis of multiple exceptional points related to three interacting eigenmodes in a 

non-Hermitian Hamiltonian 
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We have investigated the exceptional points (EPs) which are degeneracies of a non-Hermitian 
Hamiltonian, in the case that three modes are interacting with each other. Even though the para- 
metric evolution of the modes cannot be uniquely determined when encircling more than two EPs 
once, we can recover the initial configuration of the modes by encircling two EPs three times or three 
EPs twice. We confirm our expectation by numerically calculating the modes of an open quantum 
system, two dielectric microdisks, and 3x3 matrix model. 
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Topological properties of quantum mechanics such as 
Berry phase and Wilczek-Zee holonomy have drawn much 
attention in physics community [l], Q. Recently a lot 
of effort has been made to realize quantum computa- 
tion by using the topologically degenerate ground state 
of quasi-particles obeying non-Abclian statistics due to 
its robustness against external perturbations, which is 
called as topological quantum computation All 
these topological properties are based upon quantum dy- 
namics described by Hcrmitian Hamiltonians. However, 
it has been known that non-Hermitian Hamiltonians also 
exhibit non-trivial topological nature, which differs from 
that of Hermitian case in a fundamental level. 

Non-Hermitian Hamiltonians have been used to de- 
scribe open quantum systems in many areas of science. 
Since every real quantum systems are inevitably open 
otherwise no information can be acquired from it, usual 
description using the Hermitian Hamiltonian is only ap- 
proximately correct. In non-Hermitian case the eigen- 
values are complex, and the eigenstates form a non- 
orthogonal set. The complex eigenvalues have a clear 
physical meaning; the real and the imaginary part repre- 
sent the eigenenergy of a state and its decay rate, respec- 
tively. Moreover, the degeneracies of the non-Hermitian 
Hamiltonian, called as exceptional points (EPs), exhibit 
highly non-trivial characteristics compared with those of 
the Hermitian one [E|-Q- When two external parameters 
are varied so as to form a closed loop enclosing an EP 
in the parameter space, the two eigenstates, which be- 
come degenerate at the corresponding EP, are exchanged 
with each other. This is topological since it occurs only if 
the loop encloses the EP irrespective of its precise shape. 
Such a topological property of EPs has been observed 
in experiments 

d-ELO]. Recently EPs and non-Hermitian 
Hamiltonians with PT symmetry have attracted enor- 
mous interest fl~2. 13|. 

One may ask whether such a topological nature of 
EPs can be extended into more interacting modes. For 
three interacting eigenmodes, namely i, j, and k modes, 
there exist three possible combinations to form an EP 



of double degeneracies, namely EPij for the EP between 
i and j eigenmodes, EPjk for j and k, and EPki for k 
and i. Here we do not consider the triple degeneracy, 
which will be briefly discussed later. For an individual 
EP we expect the topological properties previously men- 
tioned still hold. What then happens if two or three 
EPs are encircled by one loop? Interestingly this can- 
not be uniquely determined. If both EPij and EPjk are 
encircled, there exist two possible ways to obtain the fi- 
nal result: that is, the exchange between namely i and 
j first and consequently j and k, i.e. (i,j,k) — > (k,i,j), 
or the exchange between j and k first and then i and j, 
i.e. (i, j, k) — > (j, k, i). The former differs from the latter, 
which naturally raises a question on the non-trivial struc- 
ture of multiple EPs related to three interacting eigen- 
modes. 

In this paper we show that even though the evolution of 
the eigenmodes is not uniquely determined as mentioned 
above, the eigenvalues of the holonomy matrices describ- 
ing the corresponding evolution reveal peculiar properties 
so that by encircling two (three) EPs three times (two 
times) one can recover the initial configuration of the 
eigenmodes. Our theoretical expectation is confirmed by 
numerically finding quasibound states of an open quan- 
tum system, two dielectric microdisks, and 3x3 matrix 
model. 

First we briefly explain how the non-trivial topology 
of the EP appears by using a 2x2 matrix. Consider a 
non-Hermitian matrix H described as 



H = 



E x V 
W Eo 



(1) 



where Ex, E2, V, and W are complex numbers. The 
eigenvalues are given as E± = Eq ± \fK, where Eq = 
(Ex + E 2 )/2 andA = (Ex - E 2 ) 2 /A + VW. The EP takes 
place when A = is satisfied, which implies the EP is 
codimension-2 object. Let us describe small deviation 
from the EP by a complex number in parameter space, 
namely A = 8e %< ^ , where 6(> 0) and </> are real. It is noted 
that in the Hermitian case the degeneracy occurs only if 
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both Ei = E 2 and V — W* = arc fulfilled, implying the 
degeneracy is codimcnsion-3, and the deviation A should 
be real and positive. We then continuously vary the ex- 
ternal parameters from </> to cj) + 2tt so as to encircle the 
EP in parameter space. The corresponding eigenvalues 
then rotate around the degenerate eigenvalue, E by ir 
rather than 2ir in the complex eigenvalue space due to 

Vdert -> v / (5e l ^+ 27r ) = e^VP = -VSe 1 * (2) 

implying that two eigenvalues are exchanged; {Eq + 
VE, E - VA} -> { E - VA, E + y/A} . Since eigenstates 
are uniquely determined by a given eigenvalue unless de- 
generacies occur the corresponding eigenstates have to 
be exchanged with each other as shown in Fig. 1(b). 
The so-called holonomy matrix, which is a unitary matrix 
describing the evolution of cigenmodes after parametric 
variation along a loop encircling the EP once, is expressed 
as 

M(l) = ( J J ) ■ ( 3 ) 

In order to recover the initial configuration of the eigen- 
states one should encircle the EP twice since M(l) 2 = 1. 
In fact, in addition to rearrangement of the cigenmodes 
the sign changes originating from the geometric phase 
also occurs j!,!!], for example M(l)i2 = —1 rather than 
1. However, we ignore the sign for simplicity for the time 
being. We will discuss it later. 

To analyze three interacting eigenmodes, we need to 
examine eigenvalues of 3 x 3 non-Hermitian matrix, which 
is completely described by roots of a cubic secular equa- 
tion, na mely A 3 + aX 2 + b\ + c = 0. By using Cardano's 
method [PJ] , the general solutions are given as 

{Ai,A2,A3} = {a++ev_— /3, u>a + +d>a_— /3, uja + +u>a-—/3}, 

(4) 

where a± = (g± \/ q 2 + p 3 ) , and (3 = a/3. Here 

p = 6/3 - a 2 /9, q = -c/2 + ab/6 - a 2 /27, and u = (-1 + 
V / 3«)/2. Note that w 3 = 1. Co is a complex conjugate of 
lu. The EP then takes place if a + = a-, ua + = a_, or 
w«i = a_ is satisfied, which implies the square root part 
of a± vanishes. The EPs of three interacting eigenmodes 
thus has the equivalent topological properties of those of 
two eigenmodes. It leads us to construct the holonomy 
matrices of three eigenmodes; 

/ 1 \ / 1 \ / 1 \ 

M12 =10 M 23 = 1 M 3 i =010 
\ I J \0 10/ \1 0/ 

(5) 

Make sure that [My , 'M.jk] ^ as mentioned above. 

First let us consider the case of encircling two EPs with 
a single loop. There exist six possible ways to construct 
the holonomy matrix denoted as M(2), i.e. M12M23, 
M 23 M 12 , M23M31, M 31 M 23 , M31M12, and M 12 M 31 . 




FIG. 1: (color online) Schematic diagrams describing evolu- 
tions of eigenmodes in various situations of encircling EPs. 
(a) Encircling a single EP in parameter space. Two pa- 
rameters are continuously varied so as to enclose the loop 
denoted as the black circle. The blue line attached to the 
EP represents where the avoided level crossing occurs. The 
evolution of eigenmode is described by the holonomy matrix 
M(l). (b) Two eigenmodes are exchanged by applying M(l) 
once, (c) Encircling two EPs. (d) Three eigenmodes are re- 
arranged by successive application of two different M(l)'s, 
namely My and M^i, resulting in shifting the eigenmodes, 
that is k) — > (k, i,j)- (e) Encircling three EPs. (f) Three 
eigenmodes are rearranged by applying successive application 
of three different M(l)'s, resulting in simple exchange of two 
eigenmodes, that is (i,j,k) — > (i,k,j). 



Interestingly one can obtain the unique secular equation, 
7 3 = 1, from the eigenvalue problem of all of these holon- 
omy matrices, that is, M(2)|</>) = j\4>), where \<j>) is an 
eigenstate. It tells us that if one encircles two EPs three 
times the configuration of eigenmodes returns to the ini- 
tial one since 7 3 = 1 implies M(2) 3 = 1 no matter what 
|0) is. 

As mentioned in the introduction there exist two possi- 
ble holonomy matrices describing the evolution of eigen- 
modes when encircling two EPs once for a given loop be- 
cause there are two possible orderings of two matrices. To 
decide which order is correct one should know where the 
operation performed by Mjj is completed on the path of 
the evolution in parameter space. Near the EP, roughly 
speaking, the exchange of two eigenmodes occurs when 
the avoided level crossing is passed. If we know where the 
avoided crossing takes place in parameter space, which is 
represented as a line attached to the EP in Fig. 1(a), (c) 
and (e), the holonomy matrix can be precisely expected. 



However, if the loop is far away from the EP, it is not 
easy to directly apply such an idea since the location of 
the avoided level crossing cannot be unambiguously iden- 
tified. Nevertheless, for a given loop encircling two EPs 
there exist only two holonomy matrices irrespective of its 
shape. It is noted that if the geometric phase is taken 
into account, for a given loop there exist two holonomy 
matrices depending on the direction of the rotation of 
evolution along the loop even in encircling a single EP 
as will be shown in Eq. (7) [Hj]. This also introduces 
additional non-uniqueness of the holonomy matrices. 

Next let us consider the case of encircling three 
EPs with a single loop, in which there exist 
six possible holonomy matrices denoted as M(3), 
i.e. M 3 iM 23 Mi 2 , Mi2M 3 iM 23 , M 23 Mi 2 M 3 i, 
M12M23M31, M 3 iMi 2 M 23 , and M 23 M 3 iMi 2 . Similar 
to the previous case, M(3) satisfies the secular equation, 
7 2 = 1 or 7 = 1. In order that the configuration of 
cigenmodes returns to the initial one, three EPs should 
be encircled once for a certain one dimensional solution 
space or twice for two dimensional space orthogonal to 
the previous one dimensional space. This will become 
clearer in the discussion on 3x3 matrix model later. 

We emphasize that our analysis can be directly ex- 
tended into n interacting eigenmodes, where n is an inte- 
ger (n > 3). For example, we consider four interact- 
ing modes, in which six EPs exist. If three EPs are 
enclosed, where their holonomy matrices are written as 
My, Mjfe, and M;y using the notation introduced above, 
one can show that the corresponding secular equation 
forms 7 4 = 1. It leads us to the conclusion that if three 
EPs arc enclosed four times, the initial configuration of 
the modes are recovered. 

In order to confirm our expectation we consider a phys- 
ical system described by non-Hcrmitian Hamiltonian. A 
dielectric microcavity is one of the well-known open quan- 
tum systems which have been extensively studied due to 
its wide range of application . Il7j . The nature of quasi- 
bound states or eigenmodes of electromagnetic waves 
confined in the cavities can be intuitively understood by 
considering the corresponding ray dynamics. The ray in- 
side the cavity can be totally reflected or refracted at 
the boundary depending on the angle of incidence when 
the ray hits the boundary. This clearly shows the dielec- 
tric microcavity is an open quantum system described 
by non-Hermitian Hamiltonians jlj| EH . Although there 
exist various shapes of the microcavities, here we con- 
sider two closely located microdisks, where two external 
parameters, the ratio of radii of two microdisks r and the 
distance of two microdisks d, naturally appear [2~ 
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As far as a single microdisk is concerned, the eigenmodes 
are described by two good quantum numbers, namely 
the angular quantum number m and the radial quantum 
number I. We are interested in high-Q eigenmodes whose 
corresponding ray circulates along the perimeter with a 
large angle of incidence. This is referred to as whisper- 
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FIG. 2: (color online) (a) The evolution of real eigenval- 
ues of complex wavenumbers during the parametric varia- 
tion along the closed loop enclosing two EPs, namely EPi 
and EP2 following the path A — > B — > ... — s> L — > A, 
which reveals the structure of Riemann surfaces. The closed 
loop is formed by a trapezoid with four corners, namely 
(0.750,0.260), (0.950,0.600), (0.950,0.700), and (0.750,0.360) 
in the parameter space (r, d) which corresponds to (0.2a + 
0.750,0.34a + 0.1/3 + 0.260). Two dots and the attached 
dashed lines represent EPs and the associated branch cuts, 
respectively. The level crossings and the avoided level cross- 
ing are pointed out by encircling the dashed and the solid 
circles, respectively. The avoided crossings between W2 and 
W2 modes are denoted as M and N (See the text for detail), 
(b) The shape of the eigenmodes for the eigenvalues of A to 
L designated in (a) in coordinate space. 



ing gallery eigenmode (WGM). Once two microdisks are 
considered, the WGM of one side starts to interact with 
that of the other side. However, the eigenmodes of the 
total system are dominantly localized on one of two disks 
except near the avoided level crossings 21] . It leads us 
to label the eigenmodes as WGMf^, where L(R) rep- 
resents that the eigenmode is dominantly localized in the 
left (right) disk. 
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Here we deal with four WGM's, WGM 



(6,1) 



WGM 



(7,1)' 



1.6 



WGMf 8 1): and WGMf 9 1} which are called as Wi, W 2 , 
W3, and W2 ■ The reason why we consider four rather 
than three modes is that during the parametric evolution 
W2 is exchanged with W2 , and returned since the loop 
inevitably crosses the line of avoided crossing of W2 and 
W 2 twice. Therefore, the existence of the mode W 2 can 
be ignored. We numerically find two exceptional points, 
namely EP\ related to coalescence of W\ and W3, and 
EP2 of W2 and Ws in the parameter space (r, d) . For a 
given parameter set denoted as A in Fig. 2(a), the dis- 
tributions of three eigenmodes W\, W 2 ' , and W3 in real 
space are presented in A, E and I in Fig. 2(b), respec- 
tively. Let us consider the parametric variation along a 
closed loop encircling two EPs, EP\ and EP 2 focusing 
on the evolution of W\. It is shown in Fig. 2(a) that 
the real eigenvalue of W\ evolves along the boundary of 
the surfaces denoted as the red lines during the para- 
metric variation encircling the loop three times starting 
from A. During the first encircling along the loop the 
eigenvalue follows the path A — > B — > C — > D — > E as 
shown in Fig. 2(a) transforming W\ into W 2 ' , which is 
also reproduced by applying M(2) once. Note that W2 is 
exchanged with W 2 during this evolution as mentioned 
above. However, W2 will be exchanged with W 2 again 
during the rest of the parametric evolution so that we 
can effectively ignore W2 for the whole process. Such 
an exchange between W 2 and W2 occurs at M and N in 
Fig. 2(a). If we encircle the loop once again, the eigen- 
value follows E—>F—>G—>H—>I, and W 2 is evolved 
into W3. In order to return to the initial state, W\, we 
should encircle the loop once again leading to the path 
I—>J—>K—>L—>A, that is, totally three times. This 
is exactly what we expected. 

To confirm our expectation for the case of encircling 
three EPs, we introduce the 3x3 matrix H 3 described 
by 



P 1.8, 




FIG. 3: (color online) The energy surface of the real parts of 
the eigenvalues, A, of Eq. (6) (z-axis) as varying two external 
parameters (ce,/3). Three EPs and the corresponding branch 
lines are represented as three dots and the attached dashed 
lines, respectively. 



On the other hand, if an cigenmode starts from E and 
continuously moves along the path E — > F L — > 

E, it finally returns to its initial state after encircling 
the loop twice. This again corresponds to the case of 
7 2 = 1. The similar result is obtained if one starts from 
I. We note that we can find qualitatively the same results 
obtained in Fig. 2 if the model of Eq. (6) is applied to 
the previous case of encircling two EPs. 

For simplicity the geometric phase has been ignored 
so far. Now let us briefly discuss it. It is known that 
encircling a single EP generates geometric phase so that 
Eq. ([3]) should be rewritten as 



ei — *7i S 5 
S e 2 - 172 (5 
S S e 3 - 173 



(0) 



where e\, e 2l and e 3 are functions of external parameter 
a, while 71, 72, and 73 are functions of (3. We set e\ = 
a-3, 71 = P—l, e 2 = -a+1, 72 = -/3+3, e 3 = 0, 73 — 2, 
and S = 0.4, respectively. Figure 3 shows real parts of 
three eigenvalues in the ranges of 0.4 < a < 3.5 and 
1.6 < (3 < 2.2. In the parameter space, there exist three 
EPs at (a, /3)=(1.401, 1.948), (2.072, 1.686), and (2.959, 
2.052). If an eigenmode starts from A and continuously 
evolves along the path A^B^C^D^A, it 
finally returns to its initial state after encircling the loop 
only once. This exactly corresponds to the case that the 
eigenvalue of M(3) is equal to 1, namely 7 = 1. It means 
that one eigenmode does not change at all by applying 
M(3), which is clearly shown in Fig. 1(c) by i eigenmode. 



M(l) 



-1 




1 

-1 



(7) 



depending on the rotational direction of the parametric 
variation along the loop. In this case one should encircle 
the EP four times to return to the initial configuration 
due to M(l) 2 = —1. For encircling two EPs, however, 
the geometric phases are all canceled out after cycling 
the loop three times, implying M(2) 3 = 1 that the sign 
plays no meaningful role. For encircling three EPs, one 
obtains M(3) 2 = ±1. 

Final remark is in order. Recently the chiral behaviors 
of the triple exceptional point (TEP), where three eigen- 
modes coalesce altogether, has been presented, where en- 
circling TEP shifts three eigenmodes in a cyclic manner 
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It is 



like encircling two EPs as shown above 
ascribed to the cubic-root singularity appearing in the 
triple roots of the secular equation of 3x3 matrix. At 
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the first glance, encircling three separate EPs seems to 
be similar to encircling TEP since in both cases all de- 
generacies, related to three interacting eigenmodes, are 
enclosed by a single loop. However, the eigenmode con- 
figuration returns to its initial state after going around 
degeneracies twice in the case of three separate EPs while 
three times in the case of TEP. This is not surprising be- 
cause TEP is a codimension-5 object. It is unlikely that 
the 2-dimensional surface including three separate EPs 
also contains the TEP in 5-dimensional space in general. 

In summary, we have investigated the multiple EPs re- 
lated to three interacting eigenmodes of non-Hermitian 
Hamiltonians. Encircling a certain EP of two eigenmodes 
among three has exactly the same properties of a usual 
EP of two interacting eigenmodes. When encircling more 
than two EPs once, however, the evolution of the eigen- 
modes cannot be uniquely defined. Regardless of such an 
ambiguity, the configuration of the eigenmodes return to 
their initial one when encircling two EPs three times or 
three EPs twice. We confirm our expectation for encir- 
cling two and three EPs by numerically calculating three 
interacting eigenmodes of two dielectric microdisks and 
3x3 matrix model, respectively. We believe our finding 
shed some light on study on the non-Hermitian Hamilto- 
nian. 
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